Survival of an evasive prey 

G. Oshanin * ^ O. Vasilyev P. L. Krapivsky ^, and J. Klafter ^ 

^Laboratory J.-V. Poncelet (UMI CNRS 2615), Independent University of IVIoscow, Bolshoy Vlasyevsl<iy Pereulol< 11, 119002 IVIoscow, Russia, ^ Laboratoire de PInysique 
Theorique de la Matiere Condensee (UMR CNRS 7600), Universite Pierre et Marie Curie, 4 place Jussieu, 75252 Paris Cedex 5 France,^ Max-Planck-lnstitut fur Metallforschung, 
Heisenbergstr. 3, D-70569 Stuttgart, Germany and Institut fijr Theoretische und Angewandte PInysik, University of Stuttgart, D-70569 Stuttgart, Germany,^ Department of 
Physics, Boston University, Boston, MA 02215 USA.^Sclnool of Cliemistry Tel Aviv University 69978 Tel Aviv, Israel, and H Freiburg Institute for Advanced Studies (FRIAS), 
University of Freiburg, 79104 Freiburg, Germany 

Edited by Joshua Jortner, Tel Aviv University, Tel Aviv, Israel, and approved June 24, 2009 (received for review April 20, 2009) 



We study the survival of a prey that is hunted by TV predators. The 
predators perform independent random wallcs on a square lattice with 
V sites and start a direct chase whenever the prey appears within 
their sighting range. The prey is caught when a predator jumps to 
the site occupied by the prey. We analyze the efficacy of a lazy, 
minimal-effort evasion strategy according to which the prey tries to 
avoid encounters with the predators by making a hop only when any 
of the predators appears within its sighting range; otherwise the prey 
stays still. We show that if the sighting range of such a lazy prey is 
equal to one lattice spacing, at least three predators are needed in 
order to catch the prey on a square lattice. In this situation, we es- 
tablish a simple asymptotic relation \nPcv{t) ~ (TV/V)^ In Pinjni(t) 
between the survival probabilities of an evasive and an immobile 
prey. Hence, when the density p = N/V of the predators is low, 
p <C 1, the lazy evasion strategy leads to the spectacular increase 
of the survival probability. We also argue that a short-sighting prey 
(its sighting range is smaller than the sighting range of the preda- 
tors) undergoes an effective superdiffusive motion, as a result of its 
encounters with the predators, whereas a far-sighting prey performs 
a diffusive-type motion. 

pursuit I chase | diffusion | superdiffusion [ first passage times 

Introduction and statement of results 

Pursuit-and-evasion problems have a long and fascinating 
history [1] . The classical setup involves two agents — say, 
a merchant vessel pursued by a pirate ship that it desperately 
tries to evade. The goal for both is to choose a determinis- 
tic motion strategy, given their velocities and sighting ranges, 
that optimizes their respective chances of successful pursuit 
or evasion. Similar games between adversary species occur in 
different environmental or biological systems; co-evolution of 
bacteria and phage or prey-predators contests being just two 
examples (see, e.g., [2, 3, 4] for more details). The chief differ- 
ence here is that the objects move less deterministically, their 
strategies are less "intelligent" and the number of interacting 
objects can be large. Such pursuer-evader contests are also 
assisted by some finite-range vision or smell. 



Fig. 1. A prey (a circle) and predators (crosses) on a square lattice with V sites. 



In this paper we investigate a class of pursuit-and-evasion 
problems involving a single evading prey that is being hunted 
by A'^ predators (Fig. 1) - a variation of the classic princess- 
and- monster game [5]. The predators perform independent 
nearest-neighbor random walks (RWs) on a finite square lat- 
tice (with V sites), and the prey is caught upon the first en- 
counter with a predator. In such a situation, how should the 
prey move in order to maximize its chances of not being caught 
up to time t? 

When the prey is blind, i.e. it has no information on 
predators' actual positions (its sighting range is zero), the 
best strategy is to stay still [6, 7]. The survival probability 
f'imm(i) of this immobile prey is given by 

PimmW-e-""', [1] 

where p — N/V and a depends on the diffusion coefficient of 
the predators, the structure of the lattice (particularly, the 
dimensionality), and on the number of sites V . 

The exponential decay [1] is an ultimate asymptotic for a 
finite lattice. There is also an intermediate asymptotic (which 
is the true asymptotic for infinitely large systems) [8] 

Pimm(t) =e-''^(*), [2] 

where S{t) is the mean number of distinct sites visited by a 
predator up to time t. 

The behavior of S{t) crucially depends on the spatial di- 
mension d. For nearest-neighbor RWs, S{t) grows according 
to [9] 

Si,t)-[-) and5(t)~— , [3] 

on one-dimensional (Id) and two-dimensional (2d) square lat- 
tice, respectively. In three (and higher) dimensions, S(t) 
grows linearly with time, S{t) ~ t/G, where G is the mean 
number of visits to the origin on an infinite lattice within an 
infinite time. (Hereinafter the symbol "~" signifies that we 
deal with the leading in time asymptotic behavior.) 

The question we address is how the simplest evasion strat- 
egy affects the survival probability of a prey, having a finite 
sighting range r, hunted by A'' predators with sighting ranges 
R. The predators perform independent random walks and, as 
soon as a given predator appears within distance R from the 
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prey, it changes the mode of motion and begins a direct chase, 
minimizing at each step a distance to the prey. The prey tries 
to avoid encounters with the predators by investing a minimal 
effort: Because the blind prey's best recourse is to stay still, 
we assume that the lazy prey does the same as long as all the 
predators are outside of its sighting range. Whenever preda- 
tors appear within its sighting range, the prey instantaneously 
hops to the nearest-neighboring site, chosen at random but so 
that (i) the distance from the visible predators will increase, 
and (ii) no other predator will get inside the prey's sighting 
range. We emphasize that whenever the prey hops, the choice 
of the landing site is random modulo the validity of the above 
requirements. 

We focus on the simplest situation when the sighting 
ranges of the lazy prey and of the predators are both equal to 
just one lattice spacing. We find that in this case, instead of 
obeying [l]-[2], the prey survival probability obeys: 



-Ba t 
-Bp^ S{t) 



finite lattice 
infinite lattice 



[4] 



where B is a numerical factor. The key feature is the replace- 
ment of the density p by p''. A rough explanation is that 
on the square lattice at least three predators must surround 
the prey in order to catch it. Hence, a very modest invest- 
ment in efi'ort pays back with a spectacular (several orders of 
magnitude) increase of the survival probability. 

We also observe that, due to encounters with the preda- 
tors, the prey performs long-ranged excursions on the lattice 
until the moment when it is captured. When (r,R) = (1, 1), 
the motion of the prey is effectively a diffusive motion with the 
mean-squared displacement growing linearly with time and 
the diffusion coefficients dependent on the mean density of 
predators. Surprisingly enough, there is a qualitative change 
in the behavior when R> r; we show numerically that in the 
case {r,R) = (1,2), the mean-squared displacement grows as 
t^'^^, i.e. random motion of the prey is superdiffusive. 



The model and numerical results 

The minimal model is defined as follows: 

1. There are predators on a square lattice with V sites 
and periodic boundary conditions. Predators are placed at 
random and the prey is initially at the origin of the lattice. 

2. Each predator performs a nearest-neighbor RW, the up- 
dates are made simultaneously, and the predators do not 
interact, e.g., there is no exclusion implying that a few 
predators can be at the same lattice site. 

3. If a predator is on the nearest-neighboring site to the prey, 
it hops to that site and thereby the prey is caught. 

4. Just before the predators hop, the prey checks the nearest- 
neighboring sites. If all four of them are empty, the 
lazy prey remains at the site. If some are occupied, the 
prey hops to the randomly chosen unoccupied nearest- 
neighboring site provided all four sites neighboring the tar- 
get one are unoccupied. If the latter happens, the predators 
chasing the prey still hop to the site that has been occupied 
by the prey. 

Thus, for the minimal model, the sighting ranges of both 
prey and predators are equal to one lattice spacing, r — R = 1. 
The generalization is obvious: The predator starts a direct 
chase when the distance to the prey is < R; the prey be- 
gins to move whenever there are predators within the range 
< r. By definition, the metric is Manhattan, so e.g. the site 
{x, y) is on distance \x\ + \y\ from the origin. We focus on the 



minimal model. (A few simulation results for the model with 
r = 1, R = 2 are presented below on Fig. 5.) 
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Fig. 2. Survival probability of an evasive prey on a square lattice with V = 10^ 
sites. The curves from top to bottom present the behavior of Pcv{i) when the 
density of the predators is p = 0.002, 0.003, . . . 0.009. Symbols denote results 
of numerical simulations and solid lines represent an exponential fit cxp( — b[p)t). 
Inset show/s the corresponding survival probability of an immobile prey. 



We performed Monte Carlo simulations on a 100 x 100 
square lattice with different numbers of predators. In Fig. 2, 
we plot the survival probability Pev(i) of an evasive prey and 
compare it with the survival probability Pimm(t) of an im- 
mobile prey (inset). Both display an exponential decay with 
time. The comparison reveals a pronounced difference (several 
orders of magnitude for sufficiently low density of predators) 
in the values of respective survival probabilities. In Fig. 3, 
we plot the numerical data for the characteristic decay rates 
6(p) as a function of p. This figure shows that for an immobile 
prey b{p) is a linear function of p, b(p) ~ p, as it should be, 
whereas for the lazy evasive prey b{p) ~ p^. This very dif- 
ferent scaling explains why the survival probability of evasive 
prey is so much higher when p ^ 1. 
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Fig. 3. Characteristic decay rates b(p) as a function of p. Diamonds represent 
numerical data for an immobile prey, solid line is an analytic result given by [ 1 ] with 
o = 1/G, where G is given by [10]. Circles represent numerical data for an 
evasive prey, whereas the dashed line is a fit 6(p) = 4.82 p^. 



In the following sections, we explain these findings. First, 
we consider an immobile prey and evaluate an exact expres- 
sion for Pimm(i) for finite lattices and in the thermodynamic 
limit. Then, we look on this problem from a different point of 
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view and obtain Pimm(t) by using a heuristic approach that 
reproduces correctly the temporal behavior of Pimm(t) and 
gives a correct dependence of the characteristic decay times 
on V and other pertinent parameters, but predicts slightly 
different values for numerical factors. Next, we show that in 
order for the evasive prey to be captured, the predators must 
first appear in a special configuration in which the prey will 
be under a direct attack but cannot evade. Wc specify such 
"stalemate" configurations and extend our heuristic approach 
to estimate Pev{t). 



Survival of an immobile prey 

To set up the scene, we start with a reminder on the so-called 
"target" problem [8], which describes the survival of an im- 
mobile prey in presence of predators performing independent 
RWs. This is an exactly solvable problem, but existing the- 
oretical descriptions are focused on the infinite systems. We 
determine Pimm(t) for finite V by taking advantage of some 
already known results. 

Let {Xn{t)} denote a given realization of the trajectory of 
the n**^ predator within the time interval [0,t] and X„{r) is 
its position at time r G [0,t]. We suppose here that r and t 
are integers. 

An indicator function 'P(f) of the event - that within the 
time interval [0, t] neither of the A'^ predators has visited the 
location of the prey - can be written down as 



n=l 

where P {{Xn{t)}) = 1 if X„(r) ^ for any r € [0,t] and 
otherwise. Averaging [5] we obtain^ 

Pimm{t)= [tF ^ L{X,t)\ , [6] 
\ X,X^O J 

where L {X, t) is the probability that a RW, starting at site X 
at time 0, has not yet visited the origin at time t. Explicitly, 
L{X,t) is defined (see, e.g., [10]) as 



L{X,t) = 



2ni 



dz 



t+i 1 



G(X,z) 



G{0.z) 



[7] 



where the integral is around the origin of the z plane and 
G{X, z) is the lattice Green (generating) function of the prob- 
ability G{X,t) to find a predator at site X at time t. From 
[7] we infer (see [10, 11]) that 



V - S{t), 



and consequently. 



Pimm(^) = 1 



V 



S{t) 



[81 



[91 



In the thermodynamic limit, A^ — > oo and V ^ oo with 
p = N/V kept constant, equation [9] reduces to [2]. For fi- 
nite lattices, the large time behavior of S{t) is well-known (see 
e.g., refs.[10, 12]). One finds that at sufficiently large times 



Sit) ~ V [l - e-*/^''] , G - 



2V 



ln{cV) 



d = l, 
d = 2, 



[101 



where c is a constant. For a square lattice, c « 1.8456. For 
d > 2, G is a constant as mentioned in the introduction. Over- 
all, we recover [ 1 ] with a = 1/G. 

We close this section with a remark that the target prob- 
lem is closely related to the so-called narrow escape prob- 
lem, which arises in a number of biological processes such 
as biochemical reactions in cellular microdomains (see, e.g. 
[13, 14, 15]). Here, particles (ions, molecules, protcius, etc.) 
move randomly in a bounded domain (cell, compartment) en- 
closed by a boundary that is perfectly refiecting everywhere, 
except for a small window through which particles can es- 
cape. Within this context, one is interested to calculate the 
first passage time density, which is defined as the time deriva- 
tive of the survival probability of a target placed at a position 
of the escape window in the presence of a particle diffusing in 
a bounded domain. 

In a 2d circular domain of area V with reflecting walls 
containing a small escape window of size a (size of the tar- 
get), the survival probability of such an immobile target in 
the presence of A^ particles diffusing with diffusion coefficient 
D obeys (see [14]), for V » a^, 



Pimm(t) ~ exp 



To match Pimm(t) given by [11] and our discrete-space calcu- 
lations, we set D = 1/4 and notice that the characteristic time 
of the decay in [11] is four times larger that the one defined 
by [1] with a = 1/G determined by [10], meaning that an 
immobile prey located at the origin of a periodic lattice will 
typically be caught sooner than an immobile prey hiding at 
a reflective boundary. The corresponding characteristic time 
may be even larger if the boundary is not smooth and the 
prey is located in the corner or near a cusp (see [14]). Note, 
as well, that the target size a in [11] may be interpreted here 
as the sighting range of the predators; hence, [11] signifies 
that in 2d, dependence of the characteristic decay time on the 
predators' sighting range is logarithmically weat. 



Survival of an immobile prey revisited 

We suppose from now on that f is a continuous variable, and 
predators perform continuous-time RWs with a pausing-time 
density '^(t) = exp(— t). Such a choice is for convenience only 
and would not cause any difference in the large-f behavior 
compared to discrete time RWs. 

Let Ti be a random variable defining the time within the 
interval [0,t] when the origin has been occupied by one or 
more predators simultaneously. Then, P{t) in [5] reads 



V{t) = 



1, Ti^O, 
0, Ti > 0. 



[12] 



Hence, once we know the probability density P(Ti), we may 
find the survival probability Pimm(i) of an immobile target 
from Pimm(t) = P(Ti = 0). 

Occupancy time Ti, i.e., a time spent on a given lattice 
site by one (or simultaneously k) random walkers, have been 
studied in a number of works (see, e.g. refs.[16, 17, 18]). An 
important for us result has been conjectured and verified nu- 
merically in ref.[17]; it states that P{T\) is a Gaussian 



P(Ti) ~ exp 



(Ti - Ti)' 
20-2 



[131 



^Actually the factor before the sum is equal to 1/(V — 1); since we are interested in large lattices, 
V 1, we always write V instead of V — 1. 
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and thus is entirely defined by the first two moments of Ti: 
mean value Ti and the variance = T2 - Ti . Evidently, 
Ti grows linearly with t, and hence Ti = is a large devia- 
tion from the most probable value, such that, in principle, an 
assumption that -P(Ti) is Gaussian may not be valid in this 
domain. We proceed to show, however, that an estimate 



-I imm (t) ~ P{Ti = 0) ~ exp 



[14] 



reproduces correctly the dependence of Pimm(i) on time and 
other parameters; only numerical factors in the characteristic 
decay times are incorrect. 

Let ^'o(t) be the indicator function of the event that at 
time moment r none of A'^ predators is at the origin: 



'i'o{r)=l[[l-I{X„{r))],IiX) = 



X = 
X / 0. 



Therefore, random variable Ti is explicitly defined by 



Ti = / dr (1 - *o( 
Jo 



r)). 



The first moment of Ti is now calculated to give 

/ \ N' 



Ti 



[15] 



[16] 



[17] 



where Gt ~ Gt (0) is the probability that a RW, commencing 
at the origin, is at the origin at time moment r. 
For the variance of Ti we find 



2 



dri 



dT2 <j)o,o{Ti,r2), 



[18] 



10 Jo 

where (I}o,o{ti,T2) is the two-time correlation function: 



9!>o,o(ri,T2) = 4'o(n) *o(r2) - 4'o(-ri) ■ *o(r2) 



V 



V 



V 



[19] 



Consider first the behavior predicted by [14] in the ther- 
modynamic limit. We find that 



0o,o(n,T2) ~ pG|T^_T2[, 



[20] 



implying that the occupation of the origin (or any other lat- 
tice site) of an infinite lattice is a stationary process with long- 
range algebraic correlations. Curiously enough, this behavior 
resembles very much the behavior observed experimentally for 
"blinking" of nanoscale light emitters [19]. 

Integrating [20] and noticing that Ti ~ pt, we obtain 



Ti 
2af 



fnt/\n{t), 
ft/G', 



d=l, 
d = 2, 
d > 2, 



[21] 



where /i = 37r/32, / = 1/4 and G' = drG-r. Note that G' 
has exactly the same meaning as G in [10], namely it equals 
the mean number of visits to the origin of an infinite lattice 
by a continuous-time random walk within an infiiute time. 

For finite but large V, from equations [17] and [20] we 
find that in the leading in 1/1/ order 



Ti ~ pt, ai 



2p / dri 
Jo 



dT2 G 



[22] 



which yields 




7V^/2V, 

Vln(F), 
1/G', 



d = 1, 
d = 2, 
d > 2. 



[231 



Hence, both in the thermodynamic limit and for finite V our 
approximate approach predicts correct (non trivial in low di- 
mensions) temporal evolution of Pimm (t) , correct (non-trivial) 
dependence on V, but overestimates numerical factors in the 
characteristic decay times. 



Survival of a lazy evasive prey 

We now can explain why the behavior of the survival proba- 
bility becomes so markedly difi'erent when the prey is evasive. 
First, we note the obvious fact that the evasive prey cannot 
be captured by a single predator. On a square lattice, the eva- 
sive prey cannot be captured by two predators, i.e., a more 
collective effort is required. In order to catch the prey, the 
predators have to create a stalemate type situation in which 
the prey can be attacked but cannot evade. On the square 
lattice, three predators can do this job. Several such configu- 
rations created by three predators are depicted in Fig. 4. 



X 



a 



Fig. 4. Four basic stalemate configurations from which the lazy prey cannot es- 
cape; the prey will be captured on the next step. 



When p <C 1, we may consider only lowest order "reac- 
tion events" which involve just three predators, as depicted in 
Fig. 4, and disregard higher order configurations, involving 
four and more predators. 

Due to the lazy evasion, the prey effectively performs a 
diffusive-type motion. A rigorous proof of this statement 
could be extremely challenging, and even quantifying the mo- 
tion is not trivial, as the prey is eventually caught. If we limit 
ourselves to the realizations when the lazy prey has not been 
caught up to time t and consider the mean-squared displace- 
ment of the prey, we find that it grows linearly with time (see 
Fig. 5). Its diffusion coefficient D ~ p (inset in Fig. 5) and is 
small compared with the diffusion coefficient of the predators 
(which is Dp = 1 /4, according to our definition of the motion 
of the predators). There is also another indirect argument to 
ignore the motion of the prey. In low-dimensional systems, 
the leading asymptotical behavior of the survival probability 
of the diffusive prey in presence of diffusive predators is inde- 
pendent of the diffusion coefficient of the prey [20, 21] as well 
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as of the reaction probability [11]. Hence we neglect diffusive- 
type motion of the prey and suppose that it is immobile. 
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Fig. 5. Mean-squared displacement of an evasive prey till the moment of its 
capture. Show/n are simulation results for the predator sighting ranges R = 1 
and R = 2. (Inset) The effective diffusion coefficient (determined via relation 
4lD[p) = X^[t)/t) of the prey for R. = 1. The solid line in the inset represents 
D(p) ^ p. 



Thus, the most significant effect of the lazy evasion tac- 
tics is the change of the effective reaction order — the prey 
can be captured only when any three predators appear in a 
stalemate configuration. Summing up these arguments, we 
estimate Pcv (t) of an evasive prey as 



InPev(t) ~ AC InPaW- 



[24] 



Here C = 28 is the total number of stalemate configurations, 
see Fig. 4, and P3{t) is the survival probability of an immo- 
btle prey which is captured only by three (or more) predators; 
that is, the prey can coexist with one or two predators and it 
is caught when three predators appear simultaneously at the 
site it occupies. The numerical factor A in [24] is unknown, it 
emerges since we approximately consider any stalemate config- 
uration involving three predators as a simultaneous encounter 
of predators at the same site. (One can also envision that 
stalemate configurations have different weights.) We will ex- 
tract the value of A from the fit to numerical data. 

Following our heuristic interpretation of the survival prob- 
ability of an immobile prey, we estimate Psit) as 



P3(i) ~exp 



'2a| 



[25] 



where now Ta is a random variable equal to the fraction of 
time within the interval [0, t] when the origin has been occu- 
pied by at least three predators simultaneously. 
We define next the indicator functions; 

iV iV 

*iM = E^(^'W) n [i-^(^nM)], [26] 

and 

JV JV 

*2M = E E I (Xi{r)) I {X,(r)) X 
1=1 p=i,p/i 

iV 

X n [i-nXAr))], [27] 

n— 1 ,n7^Z ,p 



of the events that at time moment r there is just one or just 
two predators at the origin, respectively. Then, T3 is given by 

Ta = / dr [1 - *o(r) - *i(r) - *2(r)] , [28] 
Jo 



and its first two moments obey 



1 N{N - l)zi / 

2 ' 



2 



[291 



and 



Tz^ = I dri f dT2 E <t>^.Ari,r2), [30] 
^0 Jo ,,^^0 



where the two-time correlation functions are defined by 

0»,,(ri,r2) = *,(ri)*,-(r2) -*I(^.17(;^. [31] 
In the leading in N /V order, Tz is simply 

7^~p't. [32] 
For finite lattices, the two-time correlation functions follow 



[331 



where ai,j{p) are polynomials in p. After straightforward, but 
tedious calculations, we find that in the leading in N/V order 
J2i j o-ijip) ~ and hence, % ~ p^J\ and cr| ^ p^ol. This 
implies that for sufficiently large times, the survival probabil- 
ity Pz{i) of an immobile prey which can be captured when any 
three predators appear simultaneously on the site it occupies, 
and the survival probability Pimm(i) of the target problem are 
related to each other through 



\nPz(t) ~ p^ lnPin,m(t). 



[341 



Consequently, for finite square lattices Pcv{t) of an evasive 
prey obeys 



InPev(t) ~ A Cp' In Pin,n.(t) ~ -AC '^^ 



ln{cV) ■ 



[351 



From the fit to numerical data, we get A « 0.55. 

In a similar fashion, we find that for an infinitely large 
square lattice 



In Pcv (t) ^-ACp^ 



ln{t)- 



[361 



Equations [35]-[36] lead to the announced result [4]. 



Discussion 

We studied the survival of a prey in the presence of predators 
performing RWs on sites of a square lattice. We analyzed a 
minimal-effort evasion tactics in which the prey tries to avoid 
predators by stepping away whenever a predator appears on 
one of the neighboring sites; otherwise the prey stays still. We 
showed that this strategy leads to the great enhancement of 
the survival probability in comparison with the stay-still strat- 
egy [6, 7]. More precisely, when the density p of the predators 
is small, the life expectancy of the immobile prey scales as 
p^^ whereas the minimal-evasion strategy results in the life 
expectancy of the order of p~^. 
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Several interesting additional conclusions are as follows: 

a) Wc assumed that the predators perform conventional 
RWs until the prey appears within their sighting ranges. 
Within this picture, the characteristic relaxation time of the 
survival probability of both an evasive and an immobile prey 
appears to be proportional to the factor G, which defines the 
mean number of returns to the origin within an infinite time 
during which a predator commences its random motion at the 
origin. In the case of a search for the prey by conventional 
RWs, in two dimensions G is large when V is large, because 
G ~ ln(y) but attains a finite value (as V —> 00) for d > 2. 
Such a behavior is associated with the fact that the spatial 
dimension d — 2 coincides with the fractal dimension of RWs. 
Playing on the side of the predators, it becomes clear that 
more efficient search for the prey in two dimensions will be 
realized when the predators' trajectories have a fractal dimen- 
sion < 2 (i.e. the best option would be to perform ballistic 
motion) . 

b) We focused on a particular situation (the minimal 
model) in which the sighting range of both the prey (r) and 
of the predators (R) were equal to one lattice spacing. When 
R > r, e.g. R — 2 and r = 1, we observed a different behav- 
ior: After the prey notices the predator, it can never escape 
from the sighting range of this predator. In this situation, as 
time increases, more and more predators turn from a random- 
search phase to the phase of the direct chase, so that the prey 
is accompanied by a tail of chasing predators. In our model, 
neither of the species is superior to the other with respect to 
speed and hence the predators who are directly chasing the 
prey are harmless - they just follow the prey but can never 
catch it, as long as the prey is not caged. We observed that in 
some realizations of the prey-predator contests on sufficiently 
small lattices, all N predators were directly chasing the prey 
which, however, survives to eternity. On the other hand, the 
predators who perform a direct chase: (i) block some possible 
directions for escape, such that it suffices now to meet just 
one randomly searching predator in order to create a stale- 
mate configuration and, (ii) exert some pressure on the prey 
prompting it to move almost ballistically. Encounters with the 



predators who have not yet started a direct chase suppresses 
pure ballistic motion, but still the prey performs effectively 
a superdiffusive motion, in contrast to the case when R = r. 
Here, the mean-square displacement X'^{t) of the prey, until 
the moment it is caught, obeys X'^{t) ~ t" with z « 1.65 
(see Fig. 5, R = 2 case) meaning that contrary to predators, 
whose random motion is recurrent in 2d and whose S{t) grows 
sublinearly with time (due to a logarithmic correction), the 
fractal dimension of the prey's trajectories ~ 1.2 < 2, and 
the mean number of distinct sites it visits grows linearly with 
t. Thus, paradoxically, although the predators who are in a 
direct chase phase are harmless to the prey, they force it to 
explore more new sites and, consequently, to die sooner. This 
interesting case bears further investigation. 

c) There should be an optimal value of the prey sighting 
range r. Clearly, a sighting range which is too small would 
not result in an effective evasion. On the other hand, having a 
too-large r is not good either. If r exceeds the mean distance 
between the predators, the prey will always see predators at- 
tempting a move and thus would most likely stay still. 

d) Finally, we remark that our results shed some light on 
an interesting problem of first passage times i to rare density 
fluctuations in diffusive systems. Within this context, one is 
interested to calculate the first passage time density Q**''(t) of 
the event when k out of A'^ particles performing RWs on a lat- 
tice appear for the first time t simultaneously on a specific lat- 
tice site [22]. Our analysis suggests that for sufficiently large 
i the first passage time density obeys lnQ^*-'(t) ~ —p*'i/G, 
with G defined by [10]. 
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